Ultimate decoherence border for matter-wave interferometry 
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Stochastic backgrounds of gravitational waves are intrinsic fluctuations of spacetime which lead 
to an unavoidable decoherence mechanism. This mechanism manifests itself as a degradation of 
the contrast of quantum interferences. It defines an ultimate decoherence border for matter-wave 
interferometry using larger and larger molecules. We give a quantitative characterization of this 
border in terms of figures involving the gravitational environment as well as the sensitivity of the 
interferometer to gravitational waves. The known level of gravitational noise determines the maximal 
size of the molecular probe for which interferences may remain observable. We discuss the relevance 
of this result in the context of ongoing progresses towards more and more sensitive matter-wave 
interferometry. 
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Fluctuations of spacetime are often referred to as a 
natural source of decoherence that defines an ultimate 
border for quantum interferences. The idea, evoked long 
ago by Feynman relies on the fact that the Planck 
mass mp = \JhcjG built on the Planck constant h, the 
velocity of light c and the Newton constant G, has a value 
~ 22 fig which lies on the borderland between microscopic 
and macroscopic masses. An object with a mass m larger 
than mp is thus associated to a Compton wavelength 
h/mc smaller than the Planck length £p = H/mpc typical 
of quantum fuzziness of spacetime. Though this length 
scale £p ~ 10 -35 m is not directly accessible to experi- 
ments, one may wonder whether fluctuation behaviours 
are modified when m crosses the mass scale mp 0, 0, Hj] • 

In this letter, we make this qualitative argument more 
specific by considering matter-wave interferometers as 
the quantum system and stochastic backgrounds of grav- 
itational waves (GW) as the source of their decoherence. 
Decoherence mechanism which might arise from Planck 
scale fluctuations of spacetime have already been stud- 
ied in the literature 0, 13 > with however a large uncer- 
tainty on the level of the latter fluctuations. Here, we 
focus our attention on known sources of spacetime fluc- 
tuations, namely GW backgrounds predicted by general 
relativity to be generated by astrophysical or cosmolog- 
ical processes. This source of noise, to be discussed 
in more detail later on, leads to an intrinsic decoher- 
ence mechanism against which interferometers cannot be 
shielded 0, IE 13 • We give a quantitative characteriza- 
tion of this mechanism in terms of relevant figures built 
up on the spectrum of the gravitational noise and the 
sensitivity of matter-wave interferometers to this noise. 

Our main purpose is to apply these ideas to the more 
and more sensitive matter- wave interferometers presently 
developed with larger and larger molecules jlfj. At 
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the moment, the decoherence of these instruments stems 
from collisions with the residual gas, emission of thermal 
radiation by the molecules or instrumental dephasings 
produced for example by vibrations of the mechanical 
structure. These noise sources can in principle be re- 
duced by using higher vacuum, lower temperature, im- 
proved velocity selection and, more generally, a better 
controlled and quieter environment available in particu- 
lar in space experiments [TlT |. With the ongoing rapid 
progress in this domain, more fundamental limits may 
eventually be reached. It is precisely the border induced 
by gravitationally induced decoherence which is investi- 
gated in the present letter. In particular, we compute 
the maximal mass for a molecular probe that preserves 
interferences. In the quantitative study presented here, 
this mass does not only depend on the Planck mass, but 
also on the geometry of the interferometer and on the 
gravitational noise level. Hence, this result brings the 
qualitative argument of Feynman to a quantitative esti- 
mation for a specific well defined physical problem. 

Besides those GW bursts which are looked for by in- 
terferometric detectors [l2j, there exist stochastic GW 
backgrounds extending over a large frequency range. A 
first part of this background is originating from the grav- 
itational emission of binary systems in our galaxy and its 
vicinity ^3- The stochastic character thus comes from 
our lack of knowledge on the precise parameters associ- 
ated with the enormous number of unresolved binaries. 
A second part of the background has a cosmological ori- 
gin coming from the primordial era of the cosmic evo- 
lution. These relic GW are produced by an amplifica- 
tion, occuring during the expansion of the Universe, of 
the primordial vacuum fluctuations of the gravitational 
field Q . For the sake of simplicity, we will suppose these 
backgrounds to be gaussian, stationary, unpolarized and 
isotropic. These simplifying assumptions are indeed suf- 
ficient for giving an estimation of decoherence induced 
by the scattering of these two backgrounds. 
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Within this context, the backgrounds are described by 
a spectral density of strain fluctuations at a given 

frequency u). This function is equivalent to the correla- 
tion of the metric in a Transverse Traceless (TT) coor- 
dinate system, at a fixed position taken here to be the 
center of the TT coordinate system : 

(M*)MQ)> = / ^ S h [u)]e- iut (1) 



Sijki = SikSn + SuS 
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Informations on the spectral density S7i[cl>] can be found 
in 0, flil UtI llij . The binary confusion background, 
deduced from the known distribution of binaries in the 
galaxy, shows a roughly flat plateau between the /iHz and 
mHz and drops rapidly on both sides on this plateau. 
The cosmological background shows al/w 3 dependence 
and should dominate at low frequencies. It depends on a 
poorly known parameter Q gw measuring the GW energy 
density compared to the critical cosmic density. 

Those GW induce a distortion of the interferomet- 
ric paths and therefore lead to a differential phase shift 
which can be evaluated simply in the eikonal approxi- 
mation [Tg| . In this approximation the Hamilton- Jacobi 
theory leads to an identification of the matter-wave phase 
<f> to the action S divided by h. The action corresponds 
to a Lagrangian density |2Jj] which couples the spatial 
part of the metric h%j to the second derivative of the 
quadrupole of the interferometer |2l| : 

Q«(t) = i J d 3 x (xV - i 5 ij x k x k J T 00 (t,x) 

This quadrupole coupling is equivalent to the dipole 
approximation used in electromagnetism to describe the 
coupling on an atom having a size much smaller than the 
wavelength. The dephasing A<£> between the two arms of 
the interferometer is given by the difference AS between 
the two action integrals and, then, by the expression 
with Q lJ replaced by the difference AQ 1 ^ between the 
quadrupoles evaluated when the probe follows either one 
arm or the other one. The resulting expression may be 
written in the frequency domain as : 



a J Zn 
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The apparatus function a tJ [oj] describes the geometry 
of the interferometer which has been assumed here to 
have a rhombic form with r the time of flight along each 
arm. Note that the previous expression of a v is valid for 
atomic interferometers where mirrors and beam splitters 



for atoms are built up on laser beams and perceived by 
atoms as freely falling objects 0- 

The expression (0 of the gravitational phase shift is 
an integral over the whole frequency spectrum. It con- 
tains two effects corresponding respectively to a global 
phase shift of the interferogram and to a reduction of its 
contrast. The separation of these two effects relies on the 
comparison of the GW frequency with the inverse of the 
measuring time T. The precise definition of this time T 
requires a complete study of a specific model of interfer- 
ometer. In the general discussion presented here, we will 
define it as the minimal time needed to build an inter- 
ferogram. T is not only larger than the time of flight r 
of the atoms along the two arms of the interferometer, 
but it has often to be much larger than t in order to 
reach a signal to noise ratio sufficient to see the interfer- 
ogram figure. While the signal to noise ratio is improved 
by averaging over a longer time T, the contrast of the 
interferogram is decreased by the change of the gravita- 
tional environment. The decoherence mechanism studied 
in this letter is precisely the result of this potential blur- 
ring of the fringes which would occur before the fringes 
have even became visible, should the noise be too large. 

Formal equivalence between the loss of contrast of in- 
terference fringes and the gen eral theory of decoherence 
has been established in |fl l22j. It turns out that the re- 
duction of the contrast exactly corresponds to the trace 
over the unobserved degrees of freedom of the gravita- 
tional environment. "Unobserved" here refers to gravi- 
tational noise lying outside the frequency window where 
fluctuations can be detected. In the following, we will de- 
note Sip the contribution of this uncontrolled noise which 
leads to a degradation of the contrast for a given signal 
processing strategy. We will also consider that this 5<p is 
deduced from A$ through a filtering function / in the 
frequency domain : 



M/>= / ^>K>]/[u 
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In the simple signal processing strategy which consists 
in averaging the interferogram over a measuring time T, 
the function f[u] is just the high pass filter defining un- 
controlled noise as frequencies larger than the inverse of 
T. More elaborated signal processing strategies could 
be studied by defining more general functions /. Deco- 
herence is then characterized by the value of the fringe 
contrast V deduced, within a gaussian description of the 
fluctuations, as the exponential of the variance of the 
uncontrolled noise |9|| : 



V = (exp(i5ip)) = exp 



A^ 2 = (V) (5) 



Using expressions ll.'il II) of the phase noises as well as 
the correlation functions Q of the metric perturbation, 
we finally rewrite the variance Aip 2 as an integral in the 
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frequency domain 



Aip 



2^ 



S h [u]A[u]F[u] 



(6) 



F[u] = |/H| 2 , ^]=Uo <i Ha B H 



The integrand is the product of three terms, the gravita- 
tional noise spectrum Sh, the apparatus response func- 
tion A and the filter function F. 

As already emphasized, A[u>] has a complicated ex- 
pression depending on the geometry of the interferom- 
eter. For forthcoming discussions, we will consider the 
commonly studied case of a Mach-Zehnder geometry 
with rhombic symmetry in the limits of small aperture 
angle (a <C 1) and non relativistic velocity (v -C c). The 
apparatus function is given by the following expression 
which captures the main ingredients of the physical de- 
scription of decoherence @ : 



-4mzH = (4f2sina) 



1 — cos(wt) 



2h 



(7) 

This response function scales as the square of the kinetic 
energy VL of the probe field measured as a frequency. We 
also remark that the function „4mz[w] goes to zero with 
the angular separation since the two arms are thus ex- 
posed to the same perturbation. Through its last term 
finally, it selects a frequency band in the gravitational 
spectrum which is essentially determined by the inverse 
of the time of flight 2r of the probe field inside the inter- 
ferometer. 

These features are sufficient to give an estimate of de- 
coherence deduced from the preceding calculations : 



Aip ~ {4flTsma) Ah 2 



(8) 



Ah 2 



dw e r l cr l ( 1 ~ cos(W) 
— S h [u)\F[u] — 
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Assuming that the measuring time T is much larger than 
r and that the filter F cuts off the potential divergence of 
the integral at its low frequency side, Ah 2 may essentially 
be interpreted as the average of the variance Ah 2 over the 
bandwidth F[u;].A[w]. 

For a first estimation, we can consider the simple as- 
sumption of a bandwidth lying within the plateau of 
the binary confusion background. This entails that Sh 
is roughly flat so that the averaged value Ah 2 is sim- 
ply given by the noise level on the plateau, that is 
St — 10~ 34 s, divided by the time of flight r. In par- 
ticular, this leads to a variance © that reproduces a 
Brownian-like diffusion of the phase characterized by a 
linear dependance of Aip 2 in the time of exposition t to 
the perturbation. In the more general discussion that 
follows, Ah 2 will be computed with the real spectrum 
of the GW backgrounds and the real bandwidth of the 
interferometer. 



We now discuss the numbers coming out of expression 
(|5|) for specific experimental configurations. Our main 
purpose is to investigate the possibility for a matter wave 
interferometer to approach the quantum/classical tran- 
sition now characterized by the quantitative condition 
Aip 2 ~ 1. Starting from the known fact that Aip 2 is usu- 
ally much smaller than unity for microscopic probes 0, 
we see on formula (JHJ that approaching Aip 2 ~ 1 requires 
two kinds of condition. Considering first the point of view 
of geometry, it is clearly needed to have an interferom- 
eter combining large angular separation and large time 
of flight. This condition is different from the large area 
condition helping the atomic interferometer to be used 
as an inertial sensor. The difference is due to the fact 
that GW vary in space and time so that even a null area 
interferometer could be sensitive to them. Moreover, sen- 
sitivity to GW is determined by the kinetic energy of the 
probe and not by its rest mass energy. This entails that 
rapid probes should be preferred to slow ones, a condi- 
tion clearly different from the one looked for with atomic 
interferometers used as inertial sensors [Tl| . 

The geometrical and energetical conditions are con- 
flicting with each other : a high velocity of the mat- 
ter beams decreases the time of flight for a given spa- 
tial size; meanwhile, it decreases the angular separation 
between interfering paths if the momentum transfer is 
fixed, which is the case for beam splitters built up on 
Raman scattering processes j^, 01 • This last argument 
can be made more precise by introducing AQ = AE^/h 
where E^ is the kinetic energy and AE^ its variation 
on the beam splitter, if we consider the case of a trans- 
fered momentum orthogonal to the velocity. With this 
notation, the variance Aip 2 is read as (4rASl) 2 Ah 2 which 
shows that the most relevant parameter for characteriz- 
ing the beam splitting is AE-^. When taking as an ex- 
ample the design of the HYPER project [2j|, AE^ has a 
very small value of the order of 10 -9 eV. This value can 
be increased by using multiple Raman scattering, up to 
140 emissions/absorptions |26j, but this is not enough for 
approaching Aip 2 ~ 1. 

The value of AE^ can be enlar ged for example by us- 
ing magnetic interaction guiding |27|. In this case, it is 
only limited by the depth of the guiding well and can 
go up to a few 10 -7 eV. Even larger values of the poten- 
tial depth (1 — 100 meV) are obtained with non resonant 
dipole interaction |2a |2jj . Deflection from material grat- 
ings also allow high kinetic energy transfer. For a slit of 
width a and a diffraction order n, the transfered momen- 
tum HAK = AE^/v has a value of the order of 2irhn/a 
which increases when a decreases. Another promising so- 
lution is the inelastic scattering of metastable molecules 
by nano-slit transmission gratings |30l | in which Van Der 
Waals interactions lead to energy transfer of the order 
AEk ~ 1 eV. For all these configurations however, the 
energy transfers are still too small to approach the tran- 
sition Aip 2 ~ 1. 
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In order to show how challenging is the objective of 
seeing the quantum/classical transition associated with 
gravitational decoherence, let us consider at this point 
an hypothetical matter-wave interferometer with a wide 
angle aperture sin a ~ 1 and, therefore, a large kinetic 
energy transfer AE^ ~ E^. This interferometer can ap- 
proach the transition Aip 2 ~ 1 if we suppose the beam 
to consist of molecules with a mass of 8 x 10 8 amu circu- 
lating at a velocity 1 km s" 1 in arms with one meter size. 
These numbers are calculated with the binary confusion 
background used as the source of fluctuations. The cos- 
mological contribution would lead to a slightly weaker 
effect if we take the value il gw — 1CP 14 . The previous 
numbers have to be contrasted with advanced projects 
of interferometers aiming at large molecules with mass 
in the 10 5 amu scale [3lJ . Meanwhile, they correspond to 
supersonic molecular beams with a kinetic energy trans- 
fer of the order of 300 keV, far above the splitting capa- 
bilities of the previously mentioned configurations. 

These numbers show that approaching the quan- 
tum/classical transition associated with gravitational de- 
coherence is out of reach for the presently developed 
molecular interferometers. An attractive idea would be 
to use Bose Einstein condensates (BEC) instead of large 
molecules. Should the BEC respond to the gravitational 
perturbation as a rigid object containing a large num- 
ber N of atoms, the parameter AE^ would have to be 
multiplied by the factor N leading to an amplification 
by a factor TV 2 of the decoherence rate. A first charac- 
terization of this rigidity condition is that the motion of 
the center of mass of the BEC induced by gravitational 
waves should correspond to frequencies lying well below 
the excitation spectrum of internal resonances. Detailed 
calculations are underway to make this characterization 
more precise. 

Anyway, this argument pleads for BEC used as an in- 
terferometric probe in space experiments where the in- 
strumental and environmental noises can be more effi- 
ciently controlled. This could be the best way to test the 
existence of an ultimate border for the observability of 
quantum interferences on matter-wave interferometers, 
due to the scattering of the spacetime fluctuations. 
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